Orbital momentum Hall effect in p-doped graphane 
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It is shown that an electric field applied to p-doped graphane generates a dissipationless orbital 
momentum Hall current. In the clean limit the corresponding Hall conductivity is independent of 
the concentration of holes. The Hall effect is related to the 27r-Berry phase accumulated when heavy 
and light holes are transported around the degeneracy point in the center of the Brillouin zone. This 
also leads to the orbital momentum edge currents in the equilibrium state, and to the accumulation 
of the orbital momentum at the edges when the system is driven out of equilibrium. 



In the last few years a great deal of attention of the 
condensed matter community has been attracted by two 
major developments: an experimental isolation of a high 
quality single- and multilayer graphene^^, and the dis- 
covery of the intrinsic spin Hall effect (SHE) in semicon- 
ductor systems with spin-orbit interaction^^. Because 
of its unique electronic properties^ graphene is consid- 
ered as one of the most prospective materials for fu- 
ture carbon-based electronics^, a possible successor of the 
standard silicon electronics. SHE is currently regarded as 
an important part of emergent spintronics^. 

Recently an important step towards graphene electron- 
ics has been made - an experimental synthesis of a fully 
hydrogenated graphene, na med gr aphane^. Graphane 
is a wide band-gap dielectri c 1 11 * 12 ! and therefore it may 
become an important part of nanoelectronic devices as 
it opens a way for a controlable creation of a potential 
landscape with spatial separation of conducting and insu- 
lating regions. However, graphane is not only a promis- 
ing material for passive parts of graphene electronics, but 
also has very interesting and exciting properties by itself. 
Theoretical works triggered by the first experiments pre- 
dict localized spin states at hydrogen vacancies^, demon- 
strate the existence of unusual strongly bound charge- 
transfer excitoniP^, and indicate that doped graphane is 
probably a high-T c superconductor!^. 

In this work I show that p-doped graphane demon- 
strates a dissipationless orbital momentum Hall effect 
(OHE). In the case of a weak disorder the correspond- 
ing Hall conductivity has a value ct oH sa e/nh indepen- 
dently of the density of holes. Similarly to the intrinsic 
SHE in systems with spin-orbit interactional, the origin 
of OHE in graphane is a Berry phase related to an effec- 
tive k-space "magnetic monopole" (Aharonov-Bohm flux 
line in this case) located at the band degeneracy point. 
In the equilibrium state it generates dissipationless or- 
bital momentum edge currents, while in the presence of 
a transport charge current it leads to the accumulation of 
the orbital momentum at the edges of the sample. These 
findings make a direct connection of graphane to the field 
of spintronics. In fact, they suggest a potential applica- 
tion of graphane as an active material for spintronics, or, 



more precisely, orbitronics^ devices, which can be natu- 
rally integrated within future graphene-based electronics 
chips. 

Graphane in its most stable chair conformatiorpJJi 2 -' is 
obtained from the ideal graphene by depositing hydro- 
gen on both sides of the graphene plane and hydrogenat- 
ing the A- and B-sublattice carbons from the upper and 
the lower side, respectively. The A (B) carbons relax 
upwards (downwards) and form cr-bonds with the corre- 
sponding H-atoms, which results in a diamond-like sp 3 
hybridization of carbon orbitals. This has a dramatic ef- 
fect on the electronic structure. Graphane is a wide, di- 
rect band-gap insulator with the gap minimum of 5.4 eV 
at the T-point (center of the Brillouin zoneJP^ 2 -'. 

States at the top of the valence band, which are 
mainly made of the carbon p-orbitals, belong to a two- 
dimensional representation E g of the graphane point 
group D 3 M Hence the corresponding k ■ p Hamiltonian 
should be a 2 x 2 matrix^. In general the symmetry fixes 
the k ■ p matrix up to a unitary transformation (choice 
of basis functions of the representation) . In Ref . [14] an 
effective Hamiltonian for holes in graphane has been de- 
rived using a set of {x 2 — y 2 , 2xy} as a basis of E g rep- 
resentation. In the present context a unitary equivalent 
form of the Hamiltonian, which corresponds to the basis 
{(x + iy) 2 , (x — ij/) 2 }, is more natural: 
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where p = — iV, p± = p x ± ip y , I is the identity 
matrix, and o± = a x ± ia y with <jj being the Pauli 
matrices. The constants 71 > 72 found irH^J deter- 
mine the effective masses of the isotropic light (L) and 
heavy (H) hole bands, m L = (7! + 72)"* ~ 0.22to and 
nT>H = (71 — 72) 1 ~ 0.64too, where mo is the bare elec- 
tronic mass. 

It is instructive to verify explicitly that H of Eq. [l] 
is indeed the proper invariant. Operations of the group 
D?,d ~ Ci x Cs v are shown on Fig. 1. These are the 
inversion i, in-plane rotations by 27r/3 about C3 axis, and 
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FIG. 1: Schematic structure of graphane plane and the sym- 
metry elements of D$d point group: the center of inversion i 
and the third order axis C3 at the center of the hexagon, and 
three vertical reflection planes a v ,i, a v ^ and cr V: s. 



reflections in three vertical planes cr v ,i, &v,2, and <j v ,s- 
Acting by the group operations on functions <fi±(x,y) = 
(x ± iy) 2 one constructs a 2x2 matrix representation of 
the group in this basiJ^ D(i) = I, 
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Using Eq. ^ we immediately find that matrices a + and 
o — belong to the representation E g and transform ac- 
cording to the rules, a± — h e ±l47r / 3 o± and a± ^""S a 



which are identical to those for the operators p± . Hence 



the term + <J-P_ in Eq. is the group invari- 

ant. Similarly the second invariant entering Eq. ([T]) is 
a product of the identity matrix I and the operator p 2 , 
which both belong to the one dimensional representation 
A\ g . The remaining linearly independent 2x2 matrix a z 
is rotationally and inversion invariant, but antisymmet- 
ric under reflections. Hence it can enter the Hamiltonian 
only in the presence of a magnetic field B, in a form of 
a z B z . It is worth noting that the Hamiltonian of Eq. ([TJ) 
is time reversal invariant with the operator of time inver- 
sion = cr x K, where K stands for the complex conjuga- 
tion. Interestingly, for 71 = Eq. JTJ coinsides with the 
Hamiltonian of bilayer graphen d 2 * lg T 

Since matrix a z acts a generator of rotations [see 
Eq. j2J|], and couples to the magnetic field as a z B z , it is 
naturally identified with the z-component of the operator 
of a local orbital momentum, i. e., the angular momen- 
tum of orbital currents inside the unit cell of graphane. 
This identification is vital for all subsequent discussion. 

Eigenstates of the Hamiltonian Eq. ([p correspond to 
two parabolic bands with dispersions E, ' = k 2 /2itil,h 
and the plane-wave spinor wave functions of the form 
f IS = $£'V kr , where 
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and k± = (k x ± ik y )/k (here k = |k[). It is customary 
to describe states of Eq. pi) as chiral 2 because they are 



eigenfunctions (with eigenvalues ±1) of a "chirality" op- 
erator <rn k , where = (cos 26* k , — sin 29k), and # k is a 
polar angle of the wave vector k. To uncover the physical 
significance of the chirality for graphane (as well as for 
bilayer graphene) we note that the chirality operator can 
be represented as follows 
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which is exactly the operator of reflection in a vertical 
plane <7„ ik parallel to the wave vector. Hence the chiral- 
ity of the L and H holes is identical to a parity, positive 
and negative, respectively, under reflection in the "prop- 
agation plane" . We will make use of this fact later on. 

An important topological property of the L and H 
states, Eq. ([3]), is a nontrivial Berry connection 
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which is related to the band degeneracy at the point 
k = 0. In fact, the degeneracy point plays a role of 
an effective Aharonov-Bohm magnetic line which carries 
a flux (Berry phase) of 2-zr 
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Now we are ready to discuss OHE in graphane. Since 
the operator of the local orbital momentum is identified 
as Lf oc = <r z we can define an operator of the orbital 
momentum current as follows 
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where V = V k -ff is the velocity operator, and {^4, B} = 
AB + BA. 

The operator of Eq. ^ is odd under reflection in the 
k-plane. Since the states of Eq. ^ are the eigenfunctions 
of this operation, all diagonal matrix elements of 3 Z are 
zero, which implies a vanishing bulk orbital momentum 
current in the equilibrium. However a nontrivial Berry 
flux suggests that a nonzero transverse current J z can 
be generated by applying an external electric field Ee iLJt . 
Using the standard linear response theory we get 



Jl = ojf (u>)E_ 
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where fg — f(E%) is the Fermi distribution function. Af- 
ter straightforward algebra (the reflection symmetry en- 
sures that only off-diagonal matrix elements contribute) 



(3) Eq. M simplifies as follows, er° H (u;) = £ijCr^ H , where 
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and A(fc) = E% - Eg = k 2 /2m L - k 2 /2m H is the en- 
ergy gap between L and H bands at a given k. Since 
for parabolic bands Ei + Eg ~ the last factor in 

Eq. ( 10 ) is proportional to the cross product of the Berry 



connection and a total derivative, Ak x Vk ln[A 2 (fc)- 



Thus at T — > the right hand side of Eq. ( 10 1 reduces 



(by partial integration) to a Fermi contour integral that 
yields, as expected, the Hall conductivity proportional to 
the Berry flux of Eq. The final result takes the form 
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where kp 



\ / 2itih£f and kp = v / 2rniiTp are the Fermi 
momenta of the H and L holes, respectively. In the limit 
uj — > the Fermi energy in Eq. (11) cancels out and we 



get the following static Hall conductivity 
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independently of the concentration of hole^H. Using the 
actual masses and recovering physical units we obtain the 
value (Tq H ~ e/irh, which closely reminds the "universal" 
intrinsic spin Hall conductivity predicted for the 2D elec- 
tron gas with spin-orbit interaction^. The crucial differ- 
ence is that the result of Eq. ( 12 ) is insensitive to a weak 
disorder - vertex corrections vanish because the Hamil- 
tonian ([lj is quadratic in the momentunPSl. The final 
static Hall conductivity in the presence of a 5-correlated 
disorder is obtained from Eq. ( 11 ) by replacing lo — > i/r, 



where r is the momentum relaxation time. 

Recently a nonzero orbital momentum Hall conductiv- 
ity has been predicted for p-doped SP^. It is instructive 
to compare OHE in graphane and in Si with intrinsic 
SHEPE While OHE in graphane is an analog of the in- 
trinsic SHE in the 2D Rashba system^ 1 , the OHE in Si 
is closely related to SHE of holes in 3D semiconductor^. 
These similarities are related to similar gauge structures 
underlying OHE and SHE in the above systems. 

A symmetry reason for the dissipationless static OHE 
is that the external electric field breaks the reflection 
symmetry in the plane perpendicular to the direction of 
E. This allows for a nonzero orbital momentum cur- 
rent propagating along that plane, i. e. OHE. Similarly 
one would expect that even in the equilibrium any re- 
flection symmetry breaking defect, e. g. the boundary, 
will generate orbital momentum currents flowing along 
the defect's potential isolines. Below I explicitly demon- 
strate the existence of such equilibrium edge currents for 
a plane boundary of graphane, modelled by an infinite 
potential wall along the line y — 0. In other words, we 
consider the "edge defect" with the potential V(r) = 
at y > 0, and V(r) = oo at y < 0. 

In the presence of an infinite, ideally reflecting poten- 
tial wall the orthonormal set of scattering states labelled 
by the incident momentum k = (k x ,—k v ), k y > takes 




FIG. 2: (Color online) Geometry of scattering and mixing of 
the light (L) and heavy (H) holes at the edge of graphane. 



the form 
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where the spinors <I>^ p and are defined by Eq. (|3|. 
The kinematics of scattering corresponding to the solu- 
tions of Eqs. ([13]) and (Thl} is illustrated on Fig. 2. The 
states v I'k( r ) an( l (rjoriginatc from incident L and 
H waves, the first terms in Eqs. (13) and (14). The sec- 



ond terms are the reflected L (H) waves with the re- 
flection coefficients _R k . Importantly, the incident L 
(H) holes also produce reflected H (L) states, the third 
terms in Eqs. (13) and (14), with amplitudes given by 

L(H) 

L— >H (H— kL) "conversion" coefficients C k . The y- 
componcnts pn (pl) of the momentum of the L— >H 
(H— kL) reflected waves are fixed by kinematic constraints 
of the conservation of energy and the a;-component of mo- 
mentum (see Fig. 2): 
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Finally, the reflection and conversion coefficients are de- 
termined from the infinite wall boundary conditions, 
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0. For the incident H-hole we find: 

x — kyPL ■ I m L 2kfr.ki, 
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(the results for the incident L-hole are obtained by inter- 
changing L and H indexes in this equation). Note that 
there is a range of incidence angles for the H-hole when 
the y-component pi, of the momentum of the reflected L- 
hole is purely imaginary, i. e., the corresponding L-state 
becomes an evanescent wave localized near the boundary. 

The L-H mixing in the scattering states Eqs. ( 13 ), ( 14 ) 
is related to the broken reflection symmetry. We can 
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see from Eq. ( 15 ) that the conversion coefficients vanish 
only for the normal incidence, k x = 0, as only in this 
case the hole states are still the eigenfunctions of the 
reflection/chirality operator, Eq. Q. 

The equilibrium orbital momentum current flowing 
along the edge is given by the expression 



(16) 
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where the prime means that the summation is restricted 
to k y > 0. Explicitly Eq. ( 16 ) can be represented in the 
following compact form 



k 



is a reduced mass, 
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From 



where 

the structure of Eq. \17\ it is clear that J^(y) decays 
away from the edge and experiences Friedel oscillations. 

For the total equilibrium edge current at T <C £f we 
get the following simple result 
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Hence the net edge current is proportional to the Fermi 
energy - the result expected on dimensional grounds as 
Ep is the only parameter of the proper dimension [h/t]. 

Now I dcmonsratc that OHE leads to the orbital mo- 
mentum accumulation at the edges of a graphane sample. 
The local orbital momentum operator a z is antisymmet- 
ric under the time reversal and the reflections. There- 
fore both symmetries must be broken to have a nonzero 
expectation value (<r z ). Let us first analyze the bulk sit- 
uation. Because of the reflection symmetry the orbital 
momentum vanishes independently for every plane wave 
state of Eq. ([3|, which in particular implies (er z ) = in 
the equilibrium. Breaking the time reversal symmetry 



and one of the reflections by applying an external elec- 
tric field E and generating a transport charge current j 
still does not produce a nonzero (<j z ). The reason is the 
remaining reflection symmetry in the vertical plane par- 
allel to j. The presence of the edge breaks that symmetry 
and thus allows for a nonzero orbital momentum, which 
should be localized at the edge as it must vanish in the 
bulk of graphane. This symmetry argumentation can be 
explicitly illustrated for the model boundary considered 
above. Because of the broken reflection symmetry the or- 
bital momentum [5 , g(r)]^cr z ^E'^(r) is nonzero for all scat- 
tering states, except for the normal incidence. However, 
it is proportional to k x (as dictated by the time rever- 
sal symmetry) and thus vanishes when integrated with 
any symmetric in k x distribution function. In the trans- 
port situation the distribution function becomes asym- 
metric and the integration yields a nonzero orbital mo- 
mentum localized at the edge. On dimensional grounds 
we conclude that the line density of the edge orbital mo- 
mentum should be proportional to H 2 j/eF ~ erE, which 
is the momentum asymmetry responsible for the trans- 
port current. This is similar to the spin accumulation at 
the surface of semiconductors described by the Luttinger 
Hamiltonian2B. 

In conclusion I predict OHE for p-doped graphane with 
all standard features of the Hall effect present: a nonzero 
dissipationlcss Hall conductivity, the equilibrium edge 
currents (which are reminiscent of the Landau diamag- 
netic currents), and accumulation of the orbital momen- 
tum at the edges. The later effect should be ovservable 
experimentally via detecting the current induced magne- 
tization of the opposite sign at the opposite edges of the 
graphane sample. 
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